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1. Introduction 

The either bosonic or fermionic character of particles available in nature fundamentally 
affects the way they behave when taking part in processes which have distinct outcomes, 
and each outcome may occur in several different manners. For example, in the celebrated 
Hong-Ou-Mandel experiment pQ, where two photons enter the input arms of a beam 
splitter (BS), the outcome when you have one photon in one output arm and one photon 
in the other arm, can occur in two different manners: Either both photons are reflected 
or are transmitted. The probability of such outcome is given by the modulus square 
of the sum of the probability amplitudes of the two occurrences. Vice versa, when the 
same experiment is performed with electrons, which are fermions as opposed to photons 
that are bosons, the outcomes are completely different: The probability of having one 
electron in one output arm and one electron in the other arm, now is given by the 
modulus square of the difference of the probability amplitudes of the two occurrences [2]. 
Also, the very quantum concept of entanglement is deeply influenced by the statistics 
of the particles [3]. An overwhelmingly part of quantum information and quantum 
communication protocols [4], quantum state tomography [5], et cetera, make use of 
multiple- way processes, entangled states and classical communication. In this respect, a 
deeper understanding of the connection between bosonic, fermionic and classical nature 
of particles is highly desirable. In the present article we study this problem and we show 
that classical and quantum statistics of identical particles are not independent concepts 
but they have a common background. 

To begin with, we consider the article "Fermion and boson beam-splitter statistics" 
[6], in which Rodney Loudon examines the action of a BS on either a pair of bosons 
or a pair of fermions. The main results of this investigation are summarized in 
table [U reproduced below. It reports: (a) the expectation values for the number 
of particles (n^) at the two output ports j E {1,2} of the BS, (b) the variance 
([An^] 2 ) = ([n^] 2 ) — (n^) 2 of the output particles number and (c) the correlation 
(n^n^) between the particles number at the two output ports. These values are 
evaluated in three different cases: for a pair of identical (i) classical]], (ii) bosonic and 
(iii) fermionic "particles"[§|, entering the BS simultaneously from opposite sides. 

From table (TJ one can see that the result reported in the second column (referring 
to identical classical particles) can be obtained as arithmetic mean of the corresponding 
results in the latter two columns (bosons and fermions). Synthetically one may 
say: identical classical particles = (bosons + fermions)/2. This result might appear 

| Think of, e. g., two white billiard balls entering a black box with two output ports. An observer who 
looks at the balls distribution at the two ports is not able to distinguish between the two balls since 
there is not any label attached to them: they are identical. However, the two particles are classical in 
the sense that if inside the black box occur different processes which yield to the same outcome A, then 
the probability of obtaining outcome A, is given by the sum of the probabilities of occurrence of each 
process leading to A. In other words, (identical) classical particles do not interfere. 
§ Being aware of the particle-wave duality, we put the word particle in quotation marks. However, for 
the sake of brevity, we will omit the quotation marks in the remainder. 
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Table 1. Expectation values (h^'), variances ([An'-'-'] 2 ) and correlation (h^fS 2 ') of 
particle numbers in the BS output arm j € {1,2}, determined for two particles that 
enter the BS simultaneously from opposite sides. 
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surprising, based on the fact that particle indistinguishability is a very different concept 
in classical and quantum physics [7] . In this work we show that this remarkable feature 
has its origin in the special structure of the two-particle Hilbert space. Furthermore, we 
generalize our observation to multi-ports (MP) (SJ |H] and to a whole class of operators. 

This article is structured as follows: in section 2 distinct two-particle states are 
introduced and their evolution is discussed. Section 3 is devoted to the description of 
identical particles. In section 4 the main result, a relation between the expectation value 
of operators invariant under particle exchange determined for a pair of identical classical 
particles, bosons and fermions, is derived. The article ends with a discussion (section 
5) and some concluding remarks (section 6). 

2. Distinct two-particle states 

For the sake of simplicity, differently from Loudon, we do not treat the problem in 
the continuous infinite Hilbert space of quantum electrodynamics, but in a finite d- 
dimensional Hilbert space. Furthermore, as generalization of the 2x2-port BS, we will 
examine an arbitrary unitary operation describing the action of a system with d input 
modes and d output modes upon a two-particle state (dxd-port, or MP for short). 

We assume that a particle located in one particular mode i G {!,... ,d} of the 
MP is described by the state vector \i). Moreover the d modes are considered to be 
fully distinguishable, i. e. (i\j) = 5ij. Hence, these states form an orthonormal basis 
of a (i-dimensional Hilbert space "H = span{|?) : i G {l,...,d}} describing all possible 
single-particle states. Here span{...} denotes the linear span of a set of vectors. [10] 

In a next step we introduce states consisting of two identical particles. We define 
these distinct two-particle states = \i) <8> \j) with \i) , \j) G % and i,j G {l,...,d}, 
where (g) denotes the tensor or Kronecker product. [H] A few words about this notation 
that is used throughout the article. The physical scenario is the following: we have 
two distinct particles, say a red (R) one and a blue (B) one, each populating one mode 
(either input or output) of a d- dimensional system. So, for example, the state where 
the red particle is occupying mode i while the blue particle is located in mode j could 
be denoted indifferently either as \R,i;B,j) or as \B,j;R,i). Then, of course, the 
kets \R,i;B,j) and \R,j;B,i) represent two physically distinct states because in the 
first case particle R is in mode i and particle B in mode j, while in the second case 
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particle R is in mode j and particle B is in mode i. When both particles occupy the 
same mode, say i, we simply could write \R,i;B,i). This somewhat heavy notation 
can be greatly lightened by choosing a specific, although arbitrary, convention for the 
ordering of the two particles. Thus, for example, from now on with \i,j) we will denote 
the state \R,i;B,j) with the convention that in the ket the first index (here i) 
always denotes the mode occupied by particle R, while the second index (here j) always 
refers to particle B. Then \i,j) and \ j,i) represent two different physical states, namely 
\R,i;B,j) and \R,j;B,i), respectively. 

Distinct two-particle states are suitable to describe any pair of fully distinguishable 
particles, such as (i) two billiard balls, (ii) a horizontally and a vertically polarized 
photon, or (iii) two electrons with opposite spin, each partner being present in mode 
% and j, respectively. Moreover, all distinct two-particle states are by definition 
orthonormal and thus do not interfere. In this sense they can be considered as describing 
classical particles. The distinct two-particle states form the standard basis of the d 2 - 
dimensional Hilbert space describing all possible two-particle states: = % <S> H = 
span{|zj) : i,j £ {1, d}}. Since has finite dimension, all operators acting on this 
Hilbert space will be bounded operators |12j . 

Under the action of the MP, an input two-particle state undergoes a linear unitary 
transformation [7 (2) : "H (2) "H (2) , defined as £/ (2) = U ' ®U with JJ :H ->H accounting 
for the linear unitary transformation of a single-particle state, i.e. \ij) = U \i)(£)U \ j) 
for all i, j G {1, d}. Fixing an ordering of the basis states by assuming \ij) to be the 
[d{i — 1) + j]th basis state, the matrix representation of is evaluated easily from the 
one of U via a Kronecker product: 



U dl U 



Ui,U 



u dd u 



(1) 



where Uij = (i\U \ j). Please note that we use the same symbol for an operator and 
its matrix representation. The probability that the MP transforms a specific input 
state |$) G %^ into a particular output state e US 2 '* is, as usual, determined by 
(tf|l7(2) |$) 2 . 

Throughout this article we will exemplify our general considerations with the help 
of two MP models. The model systems are (a) a single BS (d = 2) and (b) a combination 
of two BS (d — 3), both depicted in Fig. [TJ 

By denoting with |t| 2 the probability that a particle in input mode i is found in the 
same output mode and with |r| 2 the probability of "flipping", namely that a particle 
entering input mode i is leaving output port j, one can find that the unitary operator 
U has the following matrix representation [13] 
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Figure 1. Model systems built from BS with reflectivity r and transmissivity t. Input 
and output modes are labeled by numbers. 



Here r and t are the complex reflectivity and transmissivity coefficients of the BS. They 
do obey the constraints 



1 and rt* + tr* = 0, 



(3) 



assuring U to be unitary and thus particle number conservation. Particularly in example 
(a) a Kronecker product yields the following matrix representation of U^ 2 > = U ®U: 
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3. Symmetric and antisymmetric subspace 

The distinct two-particle states introduced in the previous section are not suitable 
for representing indistinguishable two-particle quantum states. According to quantum 
mechanics, such states are described by state vectors symmetric (bosons) or 
antisymmetric (fermions) under particle exchange. By using group theory, it is possible 
to show that the two-particle state space H^ can be expressed as direct sum of 
the symmetric subspace H s C H^ (indistinguishable bosons) and the antisymmetric 
subspace H a C (indistinguishable fermions) [15], i. e. V,^ =T-i s ® T~L a - 
Following, we define these two subspaces as: 

H s = span{|sjj) : i,j e {1, ...,d},j > i} (5a) 

= span{|s y ) : i,j E {l,...,d}}, (56) 

where \su) = \ii) and \sij) = + \ji))/y/2 for i ^ j; and 

H a = span{|aij) : i,j e {l,...,d},j > i} (6a) 

= span{|ay) : i,j G {1, ...,d},j ^ i}, (66) 

with \dij) = — \ ji))/y/2 for i ^ j. Please note that the spanning sets given in (I5~qj) 
and (16 a\) are bases of the subspaces while the ones from (156]) and (|66|) are overcomplete as 



Identical classical particles: half fermions and half bosons 



6 



\sij) = \sji) and |o^-) = — |ojj). However, we will often use the latter spanning sets since 
they simplify the notation a lot. Moreover, by definition T-L s is d(d + l)/2-dimensional 
whereas H a is d(d — l)/2-dimensional. One finds indeed =/ H s ® 1-L a and T-L s J_ H a . 
Thus, joining the two bases from (}5al) and (}6al) yields a new orthonormal basis of 7i^ 2 \ 
We will refer to it as the symmetric/antisymmetric basis. 

We can define a linear operator S : — > T-L s projecting on the symmetric part 
and a linear operator A : — > T-L a projecting on the antisymmetric part of any two- 
particle state: S = (\ij) + \ji))/2 and A\ij) = — \ji))/2 with i,j G {l,...,d}. 
Next, we give some properties of these two operators: 

• The operators S and A are projectors (S 2 = S and A 2 = A), thus possessing only 
the eigenvalues and 1. One finds that "H s (7i a ) is the eigenspace of S corresponding 
to the eigenvalue 1(0). The same property is obtained for A with the role of TL S 
and Ha swapped. 

• The operators S and A are orthogonal (SA = AS = 0) and hence do commute. We 
remark that S and A are self-adjoint operators , i.e. S = and A = A', since they 
are bounded and their matrix representation is Hermitian: (ij \ S \i'j') = S 
for all i 7 j,i',j'e{l,...,d} M- 

• The identity S + A = 1 holds and the combined operator S — A swaps the labeling 
of the two particles, i.e. 

(§-A)\ij) = \ji) with i,j G {1, d}, (7) 

as can be proved straightforwardly. 



4. Operators invariant under particle exchange 

Let us consider a linear operator O : —> that is invariant under particle 

exchange, namely an operator that satisfies the following equation: 

o = {s - Ayd(s - A), (8) 

where ([7j) should be remembered. Such operators treat the two identical particles 
forming a distinct two-particle state as indistinguishable. For example, although two 
photons identical except from orthogonal polarization are fully distinguishable, one 
cannot distinguish them with a photon detector insensitive to polarization. 
Adding (jSJ) to the trivial identity (S + Ay6(S + A) = 6, we derive 

6 = s j ds + A ] 6A. (9) 

Using this result and the properties S 2 = S, = S and the analogue for A, as well as 
SA = AS = 0, we draw two conclusions: 

• The operator O maps a state from % p always onto a state in T-L p with parity 
p G {s,a}. Consequently, O does not change the parity of symmetric and 
antisymmetric states. 

• The operator O does commute with S and A. 
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Thus, we can introduce restrictions of O to the two subspaces W s and 1-L a . The 
operator O s : rl s — > rl s shall be defined as O s = S<OS\-n s and O a : rl a — > T-L a via 
O a = A^OA\u a . Please note that the restrictions to the subspaces are necessary as the 
domain of S^OS and A^OA is whole W^. Thus, O s and O a are the equivalent of O 
in the bosonic and fermionic subspace, respectively. Additionally, the operator O, with 
help of Q), can be put into the form 

6 = &d a s + A*d a A, (io) 

since S and A are self-adjoint projectors. In summary any operator fulfilling (jBJ) consist 
of two parts acting solely either on the symmetric or antisymmetric subspace of T-L^ 2 \ 

An example of such kind of operator is the unitary transformation = U <S> U, 
because remembering that (S — A) — (S — Ay swaps the labeling of each two-particle 
state and = U <8> U, we infer that: 

(S - A)W 2 \S - A) =(S- i) f f> (2) = (S- A)\U \j) <g> U 

= U\i)®U\j) = uW\ij), (11) 

for all i,j G {1, ...,d}. Thus, due to the decomposition ( FlOl) . the matrix representation 
of in the symmetric/antisymmetric basis contains always two blocks of zeros, 
indicating forbidden transitions. The matrix representation of the restrictions Us and 
are the two non-zero blocks in the matrix representation of . In the case of our 
two model systems we find 
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for model (b), where we have defined C = y/2rt and C± = t 2 ± r 2 for the sake of 
shortness. 
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By analyzing model system (a), one sees that two identical fermions entering the 
BS simultaneously in different modes, leave it in different modes: \t 2 — r 2 \ 2 = 1 because 
of (j3J), thus obeying the Pauli exclusion principle [2]. For a 50/50 BS and two identical 
bosons in different input modes that arrive simultaneously at the BS, one obtains with 
help of (E]): \t 2 + r 2 \ 2 = 0. Therefore, they never exit the BS at different sides. This 
result is known as coalescence or Hong-Ou-Mandel effect pQ. 

Considering example (b), the matrix element (a^l ^23) vanishes since two 
identical fermions entering simultaneously mode 2 and 3 should take the same path 
when leaving the device in mode 1 and 3 but this is not permitted by Pauli's exclusion 
principle. For some transitions in model system (b) the Hong-Ou-Mandel effect can be 
observed as well, e.g., | (si 2 | \s 13 ) \ 2 = for \t\ 2 = \r\ 2 = 1/2. 

As another example of an operator that does fulfill the identity (E]), consider 

d 

n( fc ) = J2(\ik)(ik\ + \ki)(ki\), (14) 

i=l 

counting the number of particles in a particular mode k G {1, d} since: 

d 

n {k) = ^(SuSkj \lk) + 8^ \kl)) = (S kj + S ki ) ■ (15) 
1=1 

Applying h^> for d = 2, e.g., to the states 1 12) where mode 2 contains one particle and 
|22) with two particles in mode 2, yields n (2) |12) = |12) and |22) = 2 |22), as it 
should be. Note that all combinations of these number operators, like the correlation 
operator n^vP\ are invariant under particle exchange as well. A further example is 
the operator 

P^ = -^-(\kl)(kl\ + \lk)(lk\), (16) 
1 + o k i 

whose expectation value with respect to a two-particle state determines the probability 
P {ki) = (p(fcO) to find 

one particle in mode k and the other in mode I. Consider 
for example the superposition \ip) = (|11) + |12) + |21))/\/3 for d = 2. We attain 
p(u) = (^| p(n) |^) = 1/3, = p( 12 ) |^) = 2/3 and = (^| P {22) = 0. 

These two observables, namely self-adjoint operators, and p( fc ') obey a 

superselction rule since they do not connect the two subspaces % s and % a as we know 
from the general considerations above [16] . 

We obtain the following restrictions of Tv- k > to the subspaces 1-L S and H a - 

d d 

™i fc) = X^ 1 + 6ik "> \ s ik)(^k\ and n { a k) = - S ik ) \a ik )(a ik \ , (17) 

i=l i=l 

as well as of p( M ) : 

pW = \s kl )(s kl \ and P^ = Ski \aki) {a-ki \ ■ (18) 

We now have all the ingredients at hand to explain the observation stated at the 
beginning of this article. Let us consider an observable O invariant under particle 
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exchange, i.e. an operator that satisfies (jSJ), and some distinct two-particle input state 
with both particles located in different modes \ij). Subjected to the MP, the input state 
is transformed into the output density operator p = U^ 2 ^ U^ 2 \ The expectation 

value of O with respect to this output state is (O) = Tr{pO}. Recalling that the trace 
is defined as Tr{...} = J2t=i ••• l&i) f° r an arbitrary basis {\bi) : i 6 {1, d}} of H^, 
we obtain (O) = Tr s {pO} + Tr a {pO} when using the symmetric/antisymmetric basis. 
Here we denoted tracing over the symmetric and antisymmetric basis as Tr s {...} and 
Tr a {...}, respectively. In other words, we have exploited the fact = H s © "Ha to 
split the trace in into two traces performed on the subspaces T-L s and % a . Applying 
(m yields (6) = Tr s {p&6 s S} + Tr a {pA^6 a A}, where the relation S |^> = for all 
\ip) G H a and the equivalent relation for A have been used. Expressing by its 
restrictions to the subspaces H. s and 1-L a and by using orthogonality, Hermiticity and 
the projector property of S and A, we derive 

Substituting S \ij)(ij\ — \ \sij){sij\ and A A^ = ~ | we obtain one of 

the main results of this article: 

(ij\uwdu®\ij) = \ \{ Sij \ijf o s uP\ Sij ) + ( aij \u 2 Jd a u 2 a \ aij )] . (20) 

By means of this, we have proven that after passing the MP, the expectation value of 
the observable O, which is invariant under particle exchange, determined for a distinct 
two-particle state input (e.g., a pair of white billiard balls), with both partners located 
in different input modes, equals the arithmetic mean of the expectation values attained 
for the corresponding bosonic (e.g., a pair of identical photons) and fermionic (e.g., a 
pair of identical electrons) two-particle states with respect to the observables O s and 

6 a . 

If both input particles are located in the same mode, i. e. \ii), we find immediately 
U^OU^ \ii) = (su\ U 2jf d s Us 2 ^ \sa). In this case there is no difference in the 
behavior of identical bosons and identical classical particles. Note that due to Pauli's 
exclusion principle it is impossible to prepare two identical fermions simultaneously in 
the same mode. 

5. Discussion 

In table [TJ at the beginning of the article, some explicit values for the model systems 
(a) are presented. For their evaluation we considered input particles in mode 1 and 2. 
The observables are linear combinations of (for classical particles) or = ni^ 
and = hi^ (for bosonic and fermionic particles), defined in (I14p and (fTTj) . Of course 
in this case one recovers the results obtained by R. Loudon [6]. 

Applying (120]) to the operators p( kl \ pj fc ^ and from equations (fT6l) and (ITS]) , 
we obtain 

p{kl) = Ip(H) + (21) 

2 8 2 a ' 
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when the input particles are located in two different modes i and j. Thus, the transition 
probabilities at the MP between particular input and output modes for a pair of identical 
classical particles (p( fci )) equals the arithmetic mean of the transition probability for 
bosons (Pj^) and fermions (Pi ). In this spirit, we conclude that a pair of identical 
classical particles "acts" half as a pair of bosons and half as a pair of fermions. In table 
[2] we give some probabilities for model system (b) to find one particle in mode k and the 
other in mode I when the input particles where prepared in modes i and j. One verifies 
that averaging the last two columns (bosons and fermions), yields the second column 
(identical classical particles). 

Table 2. Probability to find a pair of particles in the output arms k an I, when the 
two particles entered the model system (b) simultaneously in the input modes i and j. 
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In principle this means that, for example, from determining the probability 
distribution of two unitary quantum random walks [T7J HH [19] , one performed with two 
distinguishable (e.g., orthogonally polarized) and the other with a pair of unentangled 
identical photons, the probability distribution of the same random walk for a pair of 
indistinguishable fermions can be inferred. In this spirit one can simulate the behavior of 
unentangled fermions simply with unentangled bosons. Simulating fermionic behavior 
in a quantum random walk with bosons can also be realized by using entangled states 
of polarized photon pairs \TT\ [18] . 

6. Concluding Remarks 

In summary, we have shown that the fermionic or bosonic nature of a pair of identical 
particles can be seen not only as a purely quantum feature, but it appears to be strictly 
connected to the statistical properties of classical particles. This result was achieved 
by means of a rigorous analysis of the well-established second-quantization formalism 
[20] for a two-particle, d-mode system. For more than two particles the Hilbert space 
can no longer be decomposed as direct sum of a symmetric and an antisymmetric 
subspace. Nevertheless, generalizing some of these results to more than two particles is 
certainly possible, but some more advanced group theory tools, as the Young tableaux, 
are required. 
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